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1. Introduction 

Confinement is evidently related to the nonlinearity of chromo dynamics. Since confining 
solutions of classical non-Abelian field equations are not known at present [I], it is 
believed that confinement is an essentially quantum effect. This is supported by 
numerical computations of QCD on the lattice [2j |3]. However, the analytical study 
of confinement, particularly in gauge field theory like QCD, remains a challenging task 
[3]. Thus the study of simpler field-theoretical models that simulate the characteristic 
features of confinement remains relevant. 

In this regard, it is worth mentioning earlier models, such as the dipole model [I] and 
the related higher derivative model [5] with its subsequent non-Abelian generalization 
[6]. They indicate a l/k 4 infrared behavior of the "gluon" propagator, and thus a 
linear interaction potential, even at the classical level. In spite of some quantization 
inconsistencies, these phenomenological models treat the confinement interaction as 
an elementary process, i.e., a two-particle interaction arising from the lowest-order 
approximation of perturbative dynamics of the models. 

More realistic models are the Dual Abelian Higgs model [3] and non-Abelian 
versions [3 E] in which the spontaneous symmetry breaking mechanism is used 
to generate a vacuum condensate with confining properties. In this approach the 
confinement interaction is a kind of collective effect similar to that in condensed matter 
physics. 

The two classes of models mentioned above represent quite different points of view 
on the confinement mechanism. The purpose of the present study is to shed light on 
the question: is an intermediate mechanism possible, in which confinement cannot be 
reduced to an elementary processes but is governed by cluster interactions involving 
finite numbers of particles? 

To investigate such a possibility, we utilize the variational method, in a reformulated 
Hamiltonian formalism of Quantum Field Theory (QFT), which has been demonstrated 
to be a promising and powerful approach to the relativistic bound state problem [8] 
[13]. In particular, this approach has been used to derive (and solve approximately) 
relativistic equations for two and three fermion systems, such as Positronium (Ps) and 
Muonium (Mu) [Hj, and also Ps~ and Mu" [TSj, and it was shown that the derived 
bound state energies agree with conventional perturbation theory and with experimental 
results (where available). 

The use of many-particle Fock-space components in the variational trial states leads 
to wave equations with systematically improvable bound state energy levels, as has been 
shown, for example, on the simple scalar Yukawa model [121 113] . 

In this paper, we analyze the interactions that arise from the non-linear terms in the 
mediating-field sector of the QFT Lagrangian. In particular we consider the (y9 3 +<^ 4 )- 
generalization of the Wick-Cutkosky (i.e. massless scalar Yukawa) model [16] as well as 
a version of the dipole model [H [TT] . 

We note that the models being considered are not of a non-Abelian gauge-field type. 
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The only two features which are common to the models of this paper and QCD are the 
massless and non-linear nature of the mediating field. Both features are important in 
the generating confinement but the mechanism of this effect here is different from that 
in gauge models [3l [7]. 

2. Partially reduced Wick-Cutkosky model 

The Wick-Cutkosky model [1_6] is based on the classical action integral: 

i-J*,m, (2.D 

with the Lagrangian density (h = c = 1) 

C = d^*dy - m 2 ^4> - £0*0 X + \d^X, (2-2) 

where <f>(x) is a complex scalar "matter" field with rest mass m, and X (x) is a re &l 

massless scalar field interacting with cf) via the Yukawa term g<p*<j)x (here g is an 
interaction constant). 

The stationary property of the action ( 12.ip -( !2~2l) . i.e. 51 = 0, leads to the coupled 
set of the Euler-Lagrange equations, 

(□ + m 2 )<\> = -g<f> X , (2.3) 

(□ + m 2 )0* = -^* X , (2.4) 

□X = P, (2.5) 
which determine the field dynamics; here p = —g<p*(j). 
Equation (12. 5p can be solved exactly: 

x = D*p + xo, (2.6) 
where " * " denotes the convolution [D * p] (x) = f d 4 x'D(x — x')p(x') and D(x) = 
■^5(x 2 ) is the symmetric Green function of the d'Alembert equation. If the free x 
field plays no role in the investigation the arbitrary solution xo °f t ne homogeneous 
d'Alembert equation can be omitted. Then the use of the formal solution (12. 6 j) (with 
Xo — 0) in the Lagrangian (12. 2 p leads to a self-contained variational principle for the 
interacting fields <fr(x) and <f)*(x). The modified Lagrangian £ which we shall refer to 
as the partially-reduced Lagrangian, is an important basis for the quantization of the 
model; cf. refs. [9| ITT]. 

We demonstrate here how to derive the partially-reduced Lagrangian for the Wick- 
Cutkosky model without the use of the condition xo — 0. For this purpose we consider 
the equality ( 12.61) as a change of variable x ~ > Xo where the new field xo is n °t a 'priori 
subjected to any field equation. The substitution of (12.61) directly in the Lagrangian 
(O gives 

C = 3^*8^ - m 2 0> + p(D * p + xo) + \ [dfj,{D * p + X o)d»(D * p + X o)] 
~ d (t <P*d* t <P - m 2 0> + p(D * p + xo) - \ (D * p + X o)0(D * p + X o) 
~ d^d^ - m 2 0*0 + \pD * p + ^Xo^Xo; (2.7) 
L A£f ree 
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where ~ denotes equality modulo surface terms. In this form the system is effectively 
split into two independent subsystems: the interacting <p matter field and the free Xo 
field. From this point on the physically trivial xo-dependent A£f ree term can be ignored 
(as indicated above). [j] 

The partially-reduced Lagrangian £ is non-local in space-time coordinates. The 
treatment of non-local theories of this type is a conceptually intricate, but practically 
realisable procedure. In particular, partially-reduced versions of Yukawa-like models 
are worked out in [IT] . In the next section we consider a non-linear generalization of 
Wick-Cutkosky model within the partially-reduced formulation. 



3. Nonlocal Lagrangian from a nonlinear Wick-Cutkosky model 

We proceed from the Lagrangian density 

C = d^*d^ - m 2 0> - X ~ \H<P*<P? + \d^X ~ V( X ), (3.1) 

where A > is a self-interaction coupling constant and V(x) is an arbitrary potential 
(all other quantities are the same as in (12. 2p ). 

The new terms, A(</>*0) 2 and V(x), modify the Euler-Lagrange equations (I2.3p - fl2.5p . 
In particular, the equation (I2.5P becomes the non-linear inhomogeneous d'Alembert 
equation 

□X = P-V'(X), (3-2) 

where V'(x) = dV(x)/dx- It can be formally solved by iteration (cf. ref. p2]). in the 
lst-order approximation we have: 

x = D*[p-V'(D*p)]+xo, (3.3) 

where xo includes an arbitrary solution of the homogeneous equation. 

Similarly to the case of the linear Wick-Cutkosky model, we use the replacement 
(13. 3 p (where Xo is a new field variable) in the Lagrangian (13. ip . In 1st order this gives, 

C ~ d^d^cf) - m 2 (f)*(j) +\pD*p- ±A(0*0) 2 

+ I^Xo^Xo " V(D * p + xo) + XoV'(D * p) . (3.4) 

Unlike the Lagrangian (12.71) . this functional is not completely split in the <fi and Xo 
variables. The Euler-Lagrange equation for Xo, 

D Xo = -V(D *p + X o) + V(D * p), (3.5) 

is a free-field one only in zero-order approximation. Nevertheless, it possesses the 
solution xo — which, upon substitution into (13. 4p . gives the reduced Lagrangian: 

I ~ d^d^cj) - m 2 <f)*<j) +\pD*p- iA(0*0) 2 - V(D * p) 

= C-ircc + ^i n t + ^int (3-6) 

\ It is noteworthy that, within the variational problem based on (|2.7p . the primary meaning of xo hi 
(I2.6P as general solution of the homogeneous d'Alembert equation is restored. 
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It is non-local, and the action ( 12. ip . ( 13. 6 p includes 1-, 2- and >2-fold integrations over 
Minkowsky space. 

The difference AC = C — C, i.e., the Xo-dependent part of the total Lagrangian 
( 13. 4p . is at least quadratic in the Xo variable: 

AC = AC {rcc + AC iQt where A£ free = ^d^xod^Xo, 
AC int = V(D * p) + xoV'{D *p)-V{D*p + X o) 

= ~ ^xlV"(D * p) - ^xlV"'(D * p) - . . . (3.7) 

This structure shows that the term AC is not important in the present work, as will be 
explained in more detail in Section 7. 

The non-local Lagrangian ( 13. 6 p is the lst-order approximate result of the reduction 
procedure applied to nonlinear generalizations of the Wick-Cutkosky model. In the 
Appendix we construct another local model, a kind of dipole model (with a pair of 
mediating fields), that can be reduced to the Lagrangian (I3.6P exactly. 



4. Quantization 

In order to proceed farther we need to specify the interaction potential V(x)- We choose 
V{x) = W + \*x\ (4.1) 

where k and x > are coupling constants. In this case the nonlinear Wick-Cutkosky 
model (13. ip possesses a stable perturbative vacuum and is renormalizable. 

We proceed from the partially reduced Lagrangian (I3.6p . construct the Hamiltonian 
of the model and perform the canonical quantization. Due to the non-locality of the 
Lagrangian (13. 6p . the Hamiltonization is a rather complicated procedure. It can be 
performed perturbatively, following Refs. [II] [13]. In leading-order approximation the 
Hamiltonization proceeds as follows. We work out the Hamiltonian density, 

H = Hlree + + %int + ^int ' ( 4 - 2 ) 

where 

7*!2(a0 = -2 [ dV p{x)D{x-x')p{x') + \\{<t>*{x)<t>{x)f 



1 / A^J 

2 



d x' p[x] 



D(x - x) - ^5(x 



p(x'), (4.3) 



ntl(x) = 1 K JJJ dVdV'd 4 * D(z - x)D(z - x')D(z - x")p(x)p(x')p(x"), (4.4) 

n^ t (x) = IhJJJJ dVdVdV'dSZ^ - x)D(z - x')D{z - x")D{z - x'") x 

x p(x)p{x')p(x")p(x" ! ). (4.5) 

The total interaction Hamiltonian density ( I4.2p is then expressed in terms of the Fourier 
amplitudes Ak, and A* k , B\, of the field 4>(x) (see eq. (2.14) in [13]; actually, the 
procedure is somewhat more intricate [TT] but the result is the same). Upon quantization 
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these amplitudes satisfy the standard commutation relations and become the creation 
and annihilation operators. Then the canonical Hamiltonian operator is given by 

*-/#.:«(«,.):, (4.6) 

where ": :" denotes the normal ordering of operators. Other canonical generators, such 
as linear and angular momentum, can be easily obtained. 

The term i?f ree is the standard Hamiltonian of the free complex scalar field. The 
explicit form of the pair interaction term is known (see (9j [H]) and so we shall 

(3) 

concentrate on the if; nt term. It has the following somewhat cumbersome form: 

#2 = -T^TYs [ d3fcl "- d ' fc6 J2 D(vih + V2h)D( m k 3 + m h)D( V5 k 5 + Ve k 6 ) x 
24(2vr) b J ^/k w ...kgo ^ 

"76 =± 

ii is m m vs 
x d(?7ifei + . . . + r] 6 k 6 ) :B kl Ak 2 Bk 3 Ak 4 Bk 5 Ake- , (4-7) 

where B = B, B = A\ A = A, A = B^ and the Fourier transform, D(k) = —V/k 2 , 
of the symmetric Green function of the d'Alembert equation depends on the on-shell 
4- momentum k = {/c , k}, where k = \/m 2 + k 2 . The expression (14. 7p includes 2 6 = 64 
terms. The term is of similar but more cumbersome form. We do not exhibit 
it explicitly, since, as will be seen below, it makes no contribution to the three-body 
equation derived in this work. 



5. Variational three-particle wave equations 

In the variational approach to QFT the trial state of the system is built of 
few particle channel components [T2| [13] such as the two-particle state vector 
|2) = ^/d 3 pi d 3 p 2 F 2 (pi,p 3 ) Aj^AjjO), the particle-antiparticle one |1 + 1) = 
fd 3 pi d 3 p2 G(pi,ps) A^ pi B^ 2 \0), and so on. The three-particle component has the form 

|3> = -j=J d 3 Pl d 3 p 2 d 3 p 3 Fip^p^pJAlAl^JO), (5.1) 

where the channel wave function F, which is to be determined variationally, is 
completely symmetric under the permutation of the variables P\,P2,P3- In the 
variational method the channel components, \iPn), are used to determine the 
matrix elements of the Hamiltonian, namely (ip^f\H\i/j^fi), where N,N' stand for 
1, 1, 2, 1+1, 2, 3, 2+1, 2+2,... 

(2) (3) 

We are interested here in the matrix element of the interaction H mt = + if j nt + 
of the Hamiltonian. We note that = 0, (21^12) = 0. In other 

words, purely two-particle trial states, and so the resulting variational wave equations, 

(3) 

do not sample the term H^J. Thus we first consider the three-particle case and calculate 
the matrix element 

(3|ff int |3) = /dVi-dVa d 3 p 1 ...d 3 p 3 i 7, *(p / i-P3)^(Pi-P3)/C33(p / 1 ...p^Pi...P3), (5.2) 
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where the kernel /C33 = /C 33 + /C 33 consists of the following components: 



^(Pi-Ps-Pi-Ps) = - jjT^y Kp'x + P2 + Ps -P1-P2- p 3 ) x 



X 



\M0P20P10P20 L ' 



g 2 D{p' 2 -p 2 )-\/2 , (5.3) 



,3 



^(P'l-Ps: Pl-Ps) = - 5(pi + P 2 + P 3 - Pi - P2 ~ Ps) X 

D(K-Pl)^(^ 2 -P2W 3 -P3) 



and p i0 = \J m 2 + p 2 and similarly for p'j = 1,2,3). The term H> Q ( does not 
contribute in /C 33 , i.e., K.%j = 0. 

The kernel /C33 determines the interaction in the relativistic three-particle wave 
equation that follows from the variational principle 5 (3\H — E\3) = 0, namely 

{Pw + P20 + P30 ~ E}F(p 1 ,p 2 , p 3 ) 

+ ydVidV 2 dV 3 /C33(pi,P2,P3,p / i,P 2 ,P 3 )i 71 (p , i,P 2 ,P3) = (5.5) 

where the kernel is understood to be the completely symmetrized expression (with 

respect to the variables Pi,p 2 ,p' 3 an d Pi,p 2 ,P3) of (15. 3 j) and (15. 4ft . 

(2) 

The term /C 33 of the kernel corresponds to the attractive interaction via massless 

boson exchange and repulsive contact interaction between each pair of particles while 

(3) 

/C 33 describes a cluster three-particle interaction. 

From the mathematical viewpoint the three-body wave-equation (15. 5 1) is an integral 
equation with a singular kernel. Even in simpler (say, two-particle) cases such equations 
are usually solved approximately (variationally, numerically, perturbatively), and it is 
not easy to get a general qualitative characteristic of the solutions, or to estimate the 
role of different terms of the kernel. 

In order to have some understanding of the properties of the cluster interaction we 
consider the non-relativistic limit of the equation (I5.5p . in which case the kernels simplify 
considerably, and then perform the Fourier transformation into coordinate space. In 
this representation the equation is simply a Schrodinger equation for the three-particle 
eigenfunction ^(xi, x 2 , x 3 ) (see [12]) and eigenenergy e = E — 3m: 

1 "(Pi + P\ + Pi) + V{x u x 2 , x 3 ) - el *(a?i, x 2 , x 3 ) = 0, (5.6) 



2m 

where = — iVj (i = 1, 2, 3), and the potential V(a?i, x 2 , x 3 ), like the relativistic kernel 
/C33, consists of two parts, V = V 33 + V 33 . 

v 3 y {xt,x 2 ,x 3 )= -— — - <^ -j 1 + 1 1 + 1 i 

lbirm z l\xi — x 2 \ \x 2 — x 3 \ \x 3 — X\\ 
+ {5(xi - x 2 ) + 5(x 2 - x 3 ) + 5(x 3 - cci)} , (5.7) 

V 3 f{x u x 2 ,x 3 ) = K9 U(xi,x 2 ,x 3 ). (5.8) 
(Sirm r 
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where 

U (cci, x 2 , x 3 ) = - / I n n (5.9) 

J \z - X\\\Z - x 2 \\z - x 3 \ 

The integral in r.h.s. of (15. 9p is a divergent quantity and thus equation (15. 6p may 
seem to be meaningless. However, the gradients dU(x 1 ,x 2 ,x 3 )/dx i {i = 1,2,3) which 
determine the forces in the classical background of this problem, are well defined and 
finite. Thus the "function" (15. 9 p can be presented in the form 

U (xi, x 2 , x 3 ) = U(xi, x 2 , x 3 ) + U (5.10) 

where U(x\, x 2 , x 3 ) in a regular (finite) function and Uo is an infinite negative constant 
(independent of the variables X\, x 2 , x 3 ). This constant can be absorbed by the 
eigenenergy e so that the wave equation (15. 6p gets reformulated as follows: 

3 ~ 3 2/^Q^ 

V§\xi,x 2 ,x 3 ) ^V 3 f(x 1 ,x 2 ,x 3 ) = — ~{U(x u x 2 ,x 3 ) -U } 

2/tg 3 ~ 

U(x 1 ,x 2 ,x 3 ), (5.11) 



(8nm 

e -> e = E - 3m - -^3^0 (5.12) 

where the eigenenergy e is finite (as is the potential V 33 ). 

In order to perform this reformulation explicitly, we need to resort to regularization 
of the integral (15 .9p which we consider in the next section. 

The problem of divergences is expected in the relativistic case too. But the analysis 
of the integral equation (15. 5p is a more subtle problem which shall not be undertaken 
in this work. 

6. Properties and evaluation of the 3-point potential 

Various regularization procedures are possible. In essence, one introduces some cut-off 
parameter which finally is put to (or 00). We enumerate some possibilities: 

1) We could consider the case where the mediating x field is massive, whereupon there 
would be a mass term — | /i 2 x 2 m the Lagrangian (I3.ip . In that case the gravity-like 

1 e~ Mr 

- factors would be replaced by the Yukawa forms . Thus we could regard U of eq. 

(15. 9 p as the massless-mediating-field limit of the massive-mediating-field case, 

/ e ~fi\z-xi\ e -ii\z-x 2 \ e ~n\z-x 3 \ 
d 3 z- r, (6.1) 
\Z - Xx\ \z - x 2 \ \z — x 3 \ 

which is well defined and finite for any fi > 0. 

We note that by changing the variable of integration from z to v = z — X\ in eq. 
(16. ip . we can write as 

/ e ~Liv e -^|i>+£Ei 2 | g-^|l>+a!l3| 
d 3 ^ 1 — ; ri — ; r = U^X^, X 13 ) , (6.2) 
v \v + x 12 \ « + aJi3 
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where x^ = ccj — Xj and v = \v\ 



2) Another way would be to regard U of eq. (15. 9p as a limiting R — > oo, of 



l 



E/fl(aJia, sbis) = - / dv v dv - — ■ — ■ f , (6.3) 

Jq J \V + Xi 2 \\V + X 13 \ 

where v = v/v and R is an arbitrarily large, but finite, "radius of space". 
3) We could, also, regard U of eq. (I5.9P as a limiting case, as A — > 0+, of 

U A (x 12 ,x 13 ) = - Tdv ve~ Av [ dv — 1— f . (6.4) 

Jo J \v + x 12 \\v + x 13 \ 

Evidently, any other suitable and convenient cut-off function can be used in place of 
e~ Av . 

Of course, physical results would be meaningful to the extent that they were 
independent of the choice of the regularization procedure. 

Below we establish some general properties of the regularized U function. We 
discuss in detail a convenient method of its evaluation and show that it possesses a 
logarithmic confining property when \i — > 0. 

Let us consider the regularization U fl (xi,x 2 ,x 3 ) (16. ip . It obviously obeys the 
following symmetry properties: 

(i) translational invariance: £^(aJi + A, x 2 + A, x 3 + A) = U^Xx, x 2 , x 3 ), where A £ 1R 3 ; 

(ii) rotational invariance: U^(Rxi, Rx 2 , Rx 3 ) = U^(xi,x 2 ,x 3 ), where R £ SO(3); 

(iii) permutational invariance: U^(x 2 , x\, x 3 ) = U^(xi, x 3 , x 2 ) = U^xi, x 2 , x 3 ); 

(iv) scaling transformation: U^Xxi, Xx 2 , Xx 3 ) = U\^(xi,x 2 ,x 3 ), where A £ R + . 

These properties have implications for the structure of the regularized potential. 

The properties (i)-(iii) hold for arbitrary values of the cut-off parameter /x, including 
the formal limiting case /x — > 0. Moreover, these are fundamental symmetries inherent to 
any interaction potential of a closed (nonrelativistic) system of three identical particles. 
Thus the regularized potential must possess the properties (i)-(iii) of necessity. 

The scaling property (iv) has specific implication for the regularization (I6.ip . In 
the formal limit /j — y the "function" U = U^ 1= q is scale invariant: 

(iv)' scale invariance: U(\x\, Xx 2 , Xx 3 ) = U(x\, x 2 , x 3 ), where A £ M+. 

However, as is shown below, the scaling property of the regularized potential 
U{x\, x 2 , x 3 ) is different. 

We note that an important property of the potential U(xi,x 2 ,x 3 ), with any of 
the regularizations (|6.ip - (l6.4p . follows from the symmetries 1-3, namely that it actually 
depends only on the three inter-point distances xi 2 , x± 3 , x 23 , where Xij = \xij\. Explicitly, 



e 



-/J,\v+Xij\ 



this is readily seen if the factors -j — — — — in equations (I6.ip - (l6.4p are expanded in 

spherical harmonics (fi = in (16. 3p . ( I6.4p ). the angular integrations J dv ... are carried 
out, and the orthogonality properties of the spherical harmonics are used, then (after the 
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remaining integration over dv ) , the result is seen to depend only on the lengths of the 
two vectors a3 12 , x ri and on the angle between them (or, equivalently, on x 12 , 213, £23). 

The direct calculation of the regularized potential, with any of the regularizations 
(I6.ip -( l6\4|) . is complicated. Instead, we propose a representation for the function (15. 9p 
in which its dependence on scalar arguments is manifest. This greatly simplifies the 
regularization and evaluation of U. Let us apply the well known formula: 

1 1 r 00 

- = —= I dke 

r \ 7r 



to each factor of the integrand of the expression (15. 9p (which shall be treated formally). 
Then changing the order of integration we have: 

1 



U(x 1 ,x 2 ,x 3 ) 



7T 3 /2 
■d 3 k 



k s 



_/-i.2i.2„2 1 i,2p T 2 i l,2l,2 T 2 \/j,2 
V /t l' t 2' A 12^' i 2' t 3' A 23 ^^l^i^lSI I 



dk 
~k 



(6.5) 



where k = k/k. It is obvious in this form that U(xi,X2,x 3 ) = U (x\2 , ^23^13) and, 
in addition, that the internal integral in the last line of (16. 5 j) is divergent at its lower 
boundary k — 0. 

The potential difference: 



U(x 12 ,x 2 3,xi 3 ) - U(y 12 ,y23,yn) = - dk 



dk 
~k 



e -X 2 k 2 _ e -Y 2 k 2 



(6.6) 



where X 2 = klk^x^ + k^klx^ + kfklx^, Y 2 = klk 2 yf 2 + k%kly% 3 + klk 2 y 2 3 , will be finite 
since infinite constants Uo (see (I5.10p ) from the first and second terms of (16.61) mutually 
cancel. Indeed, using the cut-off parameter e in the internal integral in r.h.s. of (16. 6 p 
yields: 



dA; 
~k 



Q -Y 2 k 2 _ Q -X 2 k 2 



00 00 



dt 

T 



Xe 

,1 1 dt 



-t 2 1 X 
— e — )-ln — , 

t e-*o Y 



Ye 



e \Ye Xe J 

i.e., the integral is convergent. 

Next, we introduce angular variables cp} on the unit sphere in fc-space, so that 
k\ = sin$cos<£>, k 2 = sin $ sin <p, k 3 — cos$. Then 

U (X12, x 23 , X13) - U (yi2, 2/23, 2/13) = W(xi2, x 23 , X13) - W(y 12 , V23, J/13), (6.7) 

where 



2% n 



W(xi2, ^13) = \ J dip J sin $ d$ In [(x 12 sin $ cos sin <£>) 2 + (a; 2 3 cos 1? sin ip) 2 

+ (x 13 cos-#cos<^) 2 ] (6.8) 
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and Xij = Xij/a . The arbitrary constant a (with dimension of length) is introduced 
so that the argument of the logarithm will be dimensionless. Actually, the potential 
difference ( 16. Tp does not depend on a while the function (16. 8 j) itself does. Since this 
function is well defined and finite, it can be considered, up to some additive constant, 
as the regularized potential: 

U (cci, x 2 , x 3 ) = W(x 12 , x 23} x 13 ) - W Q . (6.9) 

The choice of the constant Wo is a matter of taste; it can be canceled by an appropriate 
rescaling of the constant a: W(xi 2 /a, ...) = W(xi 2 /b, ...) + 47rln(6/a). Thus an 
arbitrariness of the regularized potential arises due to the scale constant a. 

We note that the regularized function (16. 9 p obeys the following scaling property: 

(iv) scale invariance: U(Xxi, Xx 2 , Xx 3 ) = U (x\, x 2 , x 3 ) + 4tt In A, where A G M+. 

The inner integration (over $) in (16. 8p can be performed explicitly. Then the change 
of variable ip — > s = cos ip yields: 

U( Xl , x 2 , x 3 ) = 4vr In Xl3 + X23 + /(£, V ), (6.10) 

4a 



where 



i 



7(4-//)= I / j. 2 arctan V / (g + e)2 + r/2 , (6.1i; 



-i 



^ _ -^13 X 23 2 _ LV^J£ 1 ■ Ay 12JL' Ay 12 ' i ■-> ' ^ ■ ' J / R ir)\ 

S — o ' " — 4 ' ^O.IZJ 

x 12 x 12 
and we have chosen for convenience: Wq = 47r(ln2 — 1). We note that the interparticle 
distances must satisfy the triangle inequalities: xi 3 + x 23 > x± 2 , x 23 + x 12 > a; 13 and 

The regularized potential ( I6.10p -( 15.12p possesses the permutational invariance (iii) 
implicitly. This is evident from the fact that any particle permutation is equivalent 
to some renumbering of k- variables in the integrals (16. 5p . ( 16 .6p and, finally, to another 
choice of angular variables in the integral ( 16. 8 p . 

In the particular cases where the points acj, x 2 and x 3 lie on a straight line the 
integral (16. lip can be calculated analytically: 

U(xx,x 2 ,x 3 ) — 4-7rln-^, where x > = max(x 12 ,Xi 3 ,x 23 ). (6.13) 

2a 

Another analytically solvable case is that of equidistant points, X\ 2 = x i3 = x 23 = r, 
whereupon in ( 16 . 12j) . £ = and rf = 3, so that I of (16. lip is a finite constant independent 
of r. Thus U(r,r,r) = 47rln(r/a) + Ci, where C\ is a finite constant, which we can ignore 
(it does not affect energy differences). For convenience we shall use "atomic units", that 

g 2 

is energies will be in units of ma 2 , and lengths in units of a = ma, where a ' 



[(xi3 + X23) 2 ~ xl 2 }[xl 2 - (xi 3 - x 23 ) 2 ] 



167T771 2 

is the dimensionless "fine structure constant" . The total potential V = V 33 + V33 (cf. 

eqs. ( 15. 7p and ( 15. 8p ) is (with A = 0), in atomic units, 

3 

V(r) = h7lnr, (r is rma, and V is V/ma 2 ) (6.14) 
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where 7 = 4k/ g. We see that V(r), in this equidistant-points subspace, is a uniformly 

3 

increasing, logarithmically confining potential (for 7 > 0). Note that V(r) ~ — for 

r 

small r (r — > + ) but V(r) ~ 7lnr for large r. Recall that if k — 7 = 0, the bound- 

3 1 

state eigenvalue spectrum (in atomic units) is the Rydberg spectrum e n = -, where 

2 n 2 

n = 1,2,3,..., and there are no bound states for e > 0. However, for 7 > 0, the 
logarithmic confining potential stretches out this Rydberg spectrum, so that there is 
a purely bound-state spectrum for e > 0. Using various approximations [T8| H~9] one 
can estimate e n ~ 7 Inn for n ^> 1. (The repulsive contact (delta-function) potentials, 
which we have ignored by taking A = 0, are of little consequence, since such repulsive 
contact potentials have an insignificant effect on the energy spectrum.) 

Other regularization methods lead, basically, to the same results. For example, 
if we use the cut-off regularization of (16 . 3 p . then for the case x 13 = 0, we obtain 
Ur = — An [1 + ln(R/x 12)] = 47rln(xi2/a) — C2 , where a is an arbitrary length parameter 
(length unit), and c 2 = 47r[l+ln(i?/a)] is a very large constant, which has to be absorbed 
into a redefined (shifted) energy, as in (|5.12p . This result is the same as eq. (I6.13p . 

In the general numerical integration of (16. lip is required. We illustrate the 

behavior of the potential in Figure 1 for the particular case X\ = a, x 2 = —a as a 
function of X3 = r. The value of potential for arbitrary configuration can be obtained 
from it using the symmetry properties (i)-(iii) and (iv). 




Figure 1. The potential U(a, —a,r) as a function of r = {x,y,z}; p = 1 x 2 + y 2 ; 
a = \a\. The function is symmetric under the inversion z — > — z and rotation around 
Oz. In particular, U = 47T0(|z| - a) In |(|z|/a + 1) if p = 0. 

In the case where one of the points is far from the others, the equality (I6.13P is 
valid asymptotically. Thus the regularized potential reveals logarithmic confinement 
properties. 

A detailed analysis of the (non-relativistic) bound-state spectrum for the general 
case requires the solution of the three body equation (15.61) . This is a quite challenging 
task in itself. However, from the confining nature of the three-point potential, we can see 
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that the spectrum will reflect confinement, much like for the equidistant-points subspace 
of fl6U4j) . 

7. Concluding remarks 

We have considered generalizations of the Wick-Cutkosky (massless scalar Yukawa) 
model that include nonlinear mediating fields. Covariant Green functions were used 
to eliminate the mediating field, thus arriving at a Lagrangian that contains nonlocal 
interaction terms. 

In the case of a massless mediating field Xi with a ^n^+^xx 4 nonlinearity, we 
evaluate the corresponding interaction term explicitly and show that the kernel has the 
form of a three- and four-point "cluster potential", cf. (I4.4p . (14. 5p . 

We consider the quantized version of this model in the Hamiltonian formalism, 
and use the variational method, with trial states built from Fock-space components, to 
derive a relativistic integral wave equation for the three-particle system. The kernels 
(relativistic potentials) are shown to contain one-quantum exchange terms and a three- 
point cluster term. In the non-relativistic limit we evaluate the explicit coordinate-space 
form of the interaction potentials and show that they consist of attractive pairwise 
Coulombic potentials and a cluster three-point confining potential. The three-point 
potential, which arises from the |kx 3 term in the Hamiltonian, is divergent (and so 
needs regularization), but the potential differences are finite. The regularized three- 
point potential is shown to be logarithmically confining, and dependent only on the 
three inter-point distances. Its evaluation, for arbitrary values of its arguments, is 
shown to be reducible to a single quadrature. 

The three-body wave-equation derived in this paper is quite complicated and must 
be solved using approximation methods. This will be the subject of forthcoming work. 

The three-particle trial state (15 .ip is found to be the simplest variational ansatz 
which manifests the confinement properties of the model. However, other sectors of 
the Fock space in the variational problem are also of interest. For example, an open 
problem is the role of the three-point interaction in the particle-antiparticle problem. 
It was pointed out in the section 5 that the simple variational particle-antiparticle trial 
state 1 1+1) does not sample the term (14.71) of the Hamiltonian. Thus, this term 
does not influence the variational wave equation derived by using only |1 + 1) (see 
[9l [HJ [12]), in which case the only Coulomb-like interaction arise. But the inclusion 
of both the |1+1) and |2+2) sectors leads to a coupled set of two many-body wave- 
equations (13] in which the effects of and are present. Whether these effects 
are confining is a question that needs to be investigated. 

Lastly, we comment on the role of "chion" Fock-space sector in the variational 
bound state problem within the reduced Hamiltonian formalism of QFT used in this 
work. This role can be examined by taking into account the Xo-dependent extra terms 
AC of the total non-local Lagrangian (13.41) . They are at least quadratic in xo including 
the free- field term ^d^Xod^Xo an d interaction terms; see eq. (13.71) . Thus the additional 
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Hamiltonian corresponding to the extra terms, AH, has no effect on variational states 
\^f) without free "chions" (i.e., quanta of the field Xo)? since (^\ AH\^f) = for such 
states. A non-trivial contribution to a variational bound-state problem may arise from 
states with two or more virtual "chions" but this is a higher-order effect in the coupling 
constants (n, x or others) of the potential V. 
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Appendix. Nonlocal Lagrangian from a nonlinear dipole model 

In this section we consider a model which is built in analogy to the linear "dipole model" 
[U [H] that simulates the confinement interaction of quarks in mesons. This model is 
nonlinear and gives Yukawa + cluster interactions. It is specified by the Lagrangian 

£ = d^d^ - m 2 0*0 - iA(0*0) 2 + p ( X + hf) + d, X - V(<p), (A.l) 

where both the xi x ) an d ^(x) are real massless scalar fields and p = —g<p*<p as in (13. ip . 

The variation of the action (I2.ip . flA.lj) leads to the coupled set of the Euler- 
Lagrange equations, 

(□ + m 2 )0 = -g 4> ( X + - X4>(4>*4>), (A.2) 

(□ + m 2 )0* = -g 0* ( X + l<p) - A0*(0», (A.3) 

Dtp = p, (A.4) 

n X = ±p-V'(<p), (A.5) 

which determine the field dynamics. 

Equations ( 1A.4I) and ( 1A.5I) possess the exact formal solution: 

(f = D*p, (A.6) 

X = D * {ip - V V)} = D * {ip - V'(D * p)} , (A.7) 

which can immediately be used in the r.h.s. of eqs. flA.2j) . flA.3j) : 

(□ + m 2 )0 = -g<f>D *{p-V\D* p)} - A0(0*0), (A.8) 

and similarly for <fi* . These equations can be derived from 51 = 0, with a Lagrangian 
identical to (13. 6p (but note that no iterative expansion, like that in eq. (13. 6p . needs to 
be made in this case). 
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